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13.1 


13.1  INTRODUCTION 

With  the  advent  of  dynamic  performance  testing  and  the  use  of  advanced 
instrumentation  and  data  reduction  techniques,  we  can  no  longer  make  the  basic 
assumptions,  such  as  wings  level  flight  or  zero  angle  of  attack,  which  allowed  the  use 
of  simplified  performance  equations.  It  is  now  necessary  to  have  a  thorough 
understanding  of  the  equations  governing  aircraft  dynamics  in  their  general  form,  and 
how  to  relate  the  information  obtained  from  instrumentation  through  these  equations. 
The  objectives  of  this  course  are  to  introduce  the  coordinate  systems  used  in  aircraft 
dynamics,  to  derive  the  coordinate  system  transformations  necessary  for  relating 
forces,  accelerations  and  velocities  acting  on  an  aircraft,  and  using  these 
transformations,  derive  the  general  velocity  and  force  relations  needed  for 
performance  data  analysis. 

13.2  AIRCRAFT  DEGREES  OF  FREEDOM  AND  THE 
EQUATIONS  OF  MOTION 

To  develop  the  equations  governing  aircraft  motion,  we  must  first  understand  what 
paths  through  space  the  aircraft  is  free  to  follow.  These  paths  are  defined  as  the 
degrees  of  freedom.  If  we  assume  the  aircraft  is  rigid,  we  find  that  in  inertial,  three- 
dimensional  space  it  has  six  degrees  of  freedom:  it  can  move  forward,  sideways,  and 
down,  and  rotate  about  any  of  the  three  orthogonal  axes.  From  Newton’s  second  law 
governing  translational  and  rotational  motion,  we  can  derive  two  vector  equations  of 
motion:  one  relating  the  external  forces  to  the  aircraft  motion  and  the  other  relating 
the  external  moments  to  the  aircraft  motion.  For  performance  testing,  we  will  only 
consider  the  force  equation.  The  moment  equation  will  be  developed  during  the 
second  half  of  this  covurse  taught  during  the  Flying  Qualities  Phase. 

Various  coordinate  systems  have  been  defined  which  allow  the  external  forces  and 
moments  to  be  easily  related.  Thus,  a  basic  understanding  of  these  coordinate 
systems  and  of  the  relationships  among  them  are  necessary  before  we  can  apply 
Newton’s  second  law  to  develop  either  of  the  equations  of  motion. 

13.3  AIRCRAFT  COORDINATE  SYSTEM  DEFINITIONS 

13.3.1  INERTIAL 

An  ineid/ial  coordinate  system  is  a  non-accelerating,  non-rotating  reference  frame  in 
which  Newton’s  second  law  is  valid.  Therefore,  the  equations  defining  the  motion  of 
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an  aircraft  must  be  developed  in  an  inertial  coordinate  system.  However,  experience 
with  physical  observations  can  be  used  to  determine  whether  a  particular  reference 
system  can  properly  be  assumed  to  be  an  inertial  coordinate  system  for  the  application 
of  Newton’s  laws  to  a  particular  problem.  For  space  dynamics  in  our  solar  system,  the 
sun  axis  system  is  a  sufficient  approximation  for  an  inertial  system.  For  aircraft 
flying  within  the  stratosphere,  the  flat  Earth  referenced  coordinate  system  is  usually 
a  sufficient  approximation  of  an  inertial  coordinate  system. 

13.3.2  FLAT  EARTH  REFERENCED 

The  flat  Earth  referenced  coordinate  system  is  a  "locally"  non-rotating  reference 
frame.  This  frame  is  made  up  of  unit  vectors  Jj,  jj,  and  which  form  an  orthogonal 
triad  defined  by  the  right-hand  rule.  The  system  is  defined  with  the  Zj-axis  ahgned 
with  the  local  gravity  vector,  toward  the  center  of  the  Earth,  leaving  the  Xi  and  Yj 
axes  to  form  the  local  horizontal  plane  (Figure  13.1).  The  Xj-axis  is  generally  aligned 
with  true  North  placing  the  Yj-axis  toward  the  East.  This  lends  the  system  to  being 
referred  to  as  the  North-East-down  coordinate  system.  The  origin  of  this  system  can 
be  fixed  at  a  particular  point  on  the  Earth’s  surface,  such  as  a  radar  site,  or  fixed  to 
the  aircraft’s  center  of  gravity.  We  will  be  interested  in  the  later  case. 


Figure  13.1  Flat  Earth  Referenced  Coordinate  System 
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13.3.3  INERTIAL  NAVIGATION  UNIT 

The  inertial  navigation  unit  (INU)  of  an  aircraft  is  designed  to  maintain  a  local 
inertial  reference  frame.  However,  the  origin  of  this  reference  frame  is  centered 
within  the  INU  which  is  not  usually  located  at  the  aircraft  center  of  gravity. 
Therefore,  during  dynamic  maneuvers,  the  INU  orientation,  acceleration  and  velocity 
data  must  be  corrected  for  fuselage  bending  and  the  offset  from  the  center  of  gravity 
in  order  to  determine  the  actual  values  in  the  true  inertial  reference  frame. 
Additionally,  the  INU  coordinate  system  may  not  be  aligned  with  the  inertial  (flat- 
Earth)  frame  defined  in  section  13.3.2.  For  example,  the  TPS  C-141A  and  F-16  INU 
provide  an  orthogonal  right-handed  inertial  system  with  the  Zj^u-axis  directed 
opposite  the  gravity  vector  (up)  and  the  Xj^u-axis  allowed  to  wander  from  true  North 
by  an  angle,  aj^u.  defined  as  the  "wander  angle"  as  shown  in  Figure  13.2.  In  contrast, 
the  F-15  aircraft  INU  provides  a  left-handed  North,  East,  up  reference  frame. 


Figure  13.2  F-16  Inertial  Navigation  Unit  Coordinate  System 
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13.3.4  BODY-FIXED 

The  body-fixed  coordinate  system  is  a  rotating  reference  frame  with  it’s  origin  fixed 
to  the  aircraft  center  of  gravity.  The  axes  form  an  orthogonal  triad  defined  by  the 
right-hand  rule  with  the  x-axis  always  pointing  through  the  aircraft’s  nose  and  the  y- 
axis  out  the  right  wing  (Figure  13.3).  This  coordinate  system  has  unit  vectors  j  g, 
and 


13.3.5  STABILITY 

The  stability  coordinate  system  is  a  rotating  reference  frame  with  the  origin  located 
at  the  aircraft  center  of  gravity.  This  system  is  made  up  of  the  unit  vectors  i^,  j^, 
and  again  forming  an  orthogonal  triad  defined  by  the  right-hand  rule.  The  x^  and 
y^-axes  form  the  plane  of  motion  of  the  aircraft,  with  the  y^-axis  always  aligned  with 
the  yg-axis  (Figure  13.4).  The  body  axes  are  related  to  the  stability  axes  through  a 
single  rotation,  defined  as  the  angle  of  attack,  a,  about  the  y^-axis.  Note  that  under 
zero  sideslip  conditions,  the  x^-axis  is  aligned  with  the  flight  path  vector  relative  to 
the  surrounding  air  mass. 
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13.3.6  WIND 

The  wind  coordinate  system  is  a  rotating  reference  frame  also  with  it’s  origin  located 
at  the  aircraft  center  of  gravity.  The  unit  vectors  2^,  j  w,  and  £y,  form  this  right- 
hand  orthogonal  system.  As  was  true  with  the  stability  coordinate  system,  the  Xyf  and 
yvv-axes  are  located  in  the  aircraft  instantaneous  plane  of  motion  (Figure  13.5). 
However,  the  x^-axis  is  always  aligned  with  the  flight  path  vector  relative  to  the 
surrounding  air  mass.  The  wind  axes  are  related  to  the  stability  axes  through  a 
single  rotation,  defined  as  the  angle  of  sideslip,  p ,  about  the  z^-axis. 


13.3.7  FLIGHT  PA-ra 

The  flight  path  coordinate  system  is  a  rotating  reference  frame  with  the  origin  located 
at  aircraft  center  of  gravity.  The  unit  vectors  2  pp,  j  pp,  and  .^pp  form  a  right-handed 
orthogonal  system.  The  x-axis  is  aUgned  with  the  flight  path  vector  relative  to  the 
inertial  reference  frame,  with  the  y-axis  perpendicular  to  the  gravity  vector  and  the 
z-axis  located  within  the  vertical  plane  formed  by  the  flight  path  and  gravity  vectors. 
The  coordinate  system  is  defined  by  two  rotations  through  the  horizontal  and  vertical 
flight  path  angles,  a  and  y,  as  shown  in  Figure  13.6.  Note  that  under  zero  wind 
conditions,  the  flight  path  vector  relative  to  the  inertial  frame  is  aligned  with  the 
flight  path  vector  relative  to  the  air  mass.  Additionally,  with  zero  winds  and  the 
aircraft  oriented  with  both  the  roll  angle,  O,  and  the  angle  of  sideshp,  p ,  at  zero,  the 
flight  path  coordinate  system  is  perfectly  aligned  with  the  wind  and  stability 
coordinate  systems.  Under  these  conditions,  y  =  0  -  a  and  a  =  T. 
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Figure  13.6  Flight  Path  Coordinate  System 


13.3.8  THRUST 

The  gross  thrust  vector,  Fg,  has  been  previously  considered  fixed  within  the  body 

coordinate  system,  at  worst  being  rotated  by  an  angle  ij-  about  the  yg-axis.  However, 
with  the  advent  of  vectored  thrust  the  gross  thrust  vector  can  now  be  rotated  in  flight 
to  produce  added  normal  and  lateral  forces  and  moments.  Therefore,  to  relate  thrust 
to  the  body-fixed  axes  now  requires  defining  a  new  coordinate  system  centered  at  the 
nozzel  exit  plane,  with  the  Xr-axis  aligned  with  the  gross  thrust  vector.  This 
coordinate  system  is  rotated  through  two  angles:  Ot  about  the  Zg-axis  and  Yt  about  an 

intermediate  y-axis,  as  shown  in  Figure  13.7.  Note  that  engine  drag,  F^,  is  assumed 
fixed  in  the  stability  coordinate  system  and  is  located  at  the  inlet  plane. 


Figure  13.7  Thrust  Coordinate  System 
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13.4  COORDINATE  SYSTEM  TRANSFORMATIONS 

With  each  of  the  aircraft  coordinate  systems  defined,  we  can  now  develop  the 
transformations  between  each  system.  This  will  allow  us  to  later  relate  force, 
moment,  or  velocity  vectors  defined  in  different  reference  frames  (Reference  1). 

13.4.1  BODY-FIXED  FROM  INERTIAL 

The  transformation  to  the  body-fixed  reference  frame  from  the  inertial  (flat  Earth) 
reference  frame  is  accomplished  by  three  rotations  through  the  angles  T,  9,  <I),  as 
shown  in  Figure  13.8.  These  angles  are  called  the  Euler  angles.  The  three  rotations 
are  not  commutative,  that  is,  they  must  be  made  in  the  specified  order  to  achieve  the 
desired  aircraft  orientation. 


X2.Xb 


Figure  13.8  Inertial  to  Body  Fixed  Coordinate  System  Transformation 

through  the  Euler  Angles 

The  first  rotation  is  about  the  Zj-axis  through  the  yaw  or  heading  angle,  T,  to  the  first 
intermediate  coordinate  system  x^,  y^,  z^,  as  shown  in  Figure  13.8a.  Thus,  to 
transform  a  vector  X  defined  in  the  inertial  reference  frame  into  coordinates  of  the 
first  intermediate  frame: 

X,  =  R,(W)  Xr 

where 

cos'?  sin'F  0 
-sin^  cos^  0 
0  0  1 


R,  ('F) 


(13.1a) 


13.8 


VOLUME  1,  PERFORMANCE  FLIGHT  TESTING 


Figure  13.8a  Inertial  to  1®*  Intermediate  Coordinate  System 
Transformation  through  the  Heading  Angle,  T 

The  second  rotation  is  about  the  yj-axis  through  the  pitch  angle,  0,  to  the  second 
intermediate  coordinate  system  Xg,  yg,  Z2,  as  shown  in  Figure  13.8b.  Defining  Yi  in 
terms  of  this  system: 

I2  =-^2(0)  II 

where 

COS0  0  -sin6 

R,(e)  =  0  1  0  (13.1b) 

sin6  0  cos6 


Figure  13.8b  1®*  to  2"^  Intermediate  Coordinate  System  Transformation 

through  the  Pitch  Angle,  0 
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The  final  third  rotation  is  about  the  Xg-axis  through  the  roll  angle,  O,  to  the  body-fixed 
coordinate  system  (Figure  13.8c).  Defining  X2  terms  of  the  body  system: 


where 


Xb  =  X2 


10  0 
0  cosO  sinO 
0  -sin^  cosO 


(13.1c) 


Figure  13.8c  2’"*  Intermediate  to  Body-Fixed  Coordinate  System 
Transformation  through  the  Roll  Angle,  O 


Now,  to  perform  the  entire  transformation  to  the  body-fixed  coordinate  system  from 
the  inertial,  we  simply  multiply  the  rotation  matrices  (equations  13.1c,  13.1b,  and 
13.1a)  to  obtain  the  desired  transformation  matrix: 


where 


X3  =  i?i(<D)  R,(Q)  R,{^)  Xr 

=  T«"Xx 


_ 


ce 

sO  s  6  cY  -  c<I>  s^P 
cOsG  cH'+s<I>s4^ 


cos'? 

sOsBs'P+cOc'P 

cd)ses'P-sd>c'F 


-s6 

s4>c6 

c<I>c6 


(13.2a) 


Note:  c  =  cos  and  s  =  sin 


J 
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Because  they  are  orthogonal,  each  of  the  matrices  defined  in  the  above  equations  have 
the  unique  property  that  the  transpose  of  the  matrix  is  equal  to  the  inverse. 

[T]T  =  [JJ-l 

Thus,  if  we  wish  to  transform  a  vector  expressed  in  body-fixed  coordinates  into  the 
inertial  frame: 

Ix  = 

where 

c6c'F  sOsOcT-  cO  s'?  c<I)  s  0  +  s3> 

=  cos'?  sO  s  0  s'F  +  c^>  c^F  c<^  s  0  s'F  -  sO  c'F  (13.2b) 

-s0  s^c0  cOc0 


13.4.2  BODY-FEOED  FROM  STABILITY 

The  transformation  to  the  body-fixed  reference  frame  from  the  stability  includes  a 
single  rotation  about  the  ys-axis  through  the  angle  of  attack,  a,  as  shown  in  Figure 
13.9.  To  transform  X  expressed  in  the  stability  coordinates  into  body  coordinates: 


X.  =  T»“  J, 

where 

cosa  0  -sina 

=  1^2  (a)  =  0  1  0  (13.3a) 

sina  0  cosa 


Figure  13.9  Stability  to  Body-Fixed  Coordinate  System  Transformation 

through  the  Angle  of  Attack,  a 
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To  transform  X  defined  in  body-fixed  coordinates  to  stability  coordinates: 


where 


Is  = 


T^®  =  R^i-a) 


cosa  0  sina 
0  10 
-sina  0  cosa 


(13.3b) 


13.4.3  WIND  FROM  STABILITY 

The  transformation  to  the  wind  reference  frame  from  the  stability  includes  a  single 
rotation  through  the  angle  of  sideslip,  p ,  as  shown  in  Figure  13.10.  To  transform^ 
expressed  in  the  stability  coordinates  into  wind  coordinates: 


Figure  13.10  Stability  to  Wind  Coordinate  System  Transformation 
through  the  Angle  of  Sideslip,  p 
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To  transform  X  defined  in  wind  coordinates  to  stability  coordinates: 


and 


Js  = 


cos  P 
sinP 
0 


-sinP  0 
cos  p  0 
0  1 


(13.4b) 


13.4.4  WIND  FROM  BODY-FIXED 

The  transformation  to  wind  coordinates  from  the  body-fixed  system  includes  the  two 
rotations  through  a  and  p  defined  in  paragraphs  13.4.2  and  13.4.3,  as  shown 
previously  in  Figure  13.5. 

To  transform  X  expressed  in  body  coordinates  into  wind  coordinates: 


Jy  =  Xb  =  Xb  =  -R3  (P  )  ^2  ( -a)  Xb 


where 


cosacosp  sinP  sinacosP 

-cosasinP  cosP  -sinasinP 

-sina  0  cosa 


To  transform  X  back  into  body  coordinates: 


X_  'J'BW  -Y"  =s  T  Y" 

B  A/W  la  J  Apv 


where 


J'BW  _ 


cosa  cos  P 
sinP 

[sina  cos  P 


-cosa  sinp 
cos  P 

-sina  sinP 


-sina 

0 

cosa 


(13.5a) 


(13.5b) 
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13.4.5  WIND  FROM  INERTIAL 

With  the  development  of  aU  transformations  between  each  coordinate  system 
complete,  the  most  general  form  of  the  transformation  matrices  between  the  inertial 
coordinate  system  and  the  wind  coordinate  system  can  now  be  built.  Using  equations 
13.2a  and  13.5a,  the  transformation  to  wind  from  inertial  is  found  to  be: 

I  =:  yBJ 


c  0  cTca  cp 

c0  sT ca  cP 

-s0  ca  cp 

+  ( s<I>  s  0  cT  -  c<I>  s'?)  s 

P  +  (sOs0  s'F+ cOc'}'}  sP 

+sO  c0  sP 

+  ( c<[>s0  o'?  +  s<I>s'F)  sa 

cP  +  ( c3>s6  s'F  -  sO  cT)  sa  cP 

+cO  c0  sa  cp 

-c0  o'?  ca  sP 

-c0  sT  ( ca sP) 

s0  ca  sp 

+  ( sO  s  0  cT  -  c<I>  s'F)  c 

P  +(s<I>s0  s'F+ c<I)cT)  cp 

+s3) c0  cP 

-  ( c<E> s 0  c'P  +  sO s'F)  sa 

sp  -  ( c<I>s0  s^' -  s^  cT)  sa  sp 

-c<I> c0  sa sP 

-c0  c'Fsa 

-c0  s4^sa 

s0  sa 

L  +  ( c<I>  s  0  s^)  ca  +  ( cQ>  s  0  -  s3>  c'F)  ca 

+c<I>  c0  ca 

(13.6a) 

And  the  transformation  to  the  inertial  from  the  wind  system: 

c0  c'F  ca  cp 

-c0  s'?  ca  sP 

-c0  c'Fsa 

+  ( s<I>  s  0  c'F  -  cO  s'?)  s  P 

+  (sa>s0  c'P- cOs'P)  cp  +(ca>s0  cT+ SOS'?) 

( c<I> s  0  o'?  +  sO s'F)  sa  cP 

-  ( c<I>  s  0  c'P  +  sO  s'?)  sa  s P 

c0  s'F  ca  cp 

-c0  s'?  (casP) 

-c0  sTsa 

+  ( s  0  sT  +  c<l>  c'F)  s  P 

+  (sa)s0  sT+ c^c'P)  cp  +(c<I)s0  s'?- sOc^) 

( c<E> s 0  sT  -  s<I)  cT)  sa  cP 

-  ( c^>  s  0  s'?  -  sO  c'P)  sa  s  P 

-s0  ca  cp 

s0  ca  sP 

s0  sa 

+s<I>c0  sP 

+s3>c0  cP 

+c<I>c0  ca 

+cOc0sacP 

-c<I>c0  sasP 

(13.6b) 


13.14 


VOLUME  1,  PERFORMANCE  FLIGHT  TESTING 


13.4.6  FLIGHT  PATH  FROM  INERTIAL 


The  transformation  from  inertial  reference  frame  to  the  flight  path  reference  frame 
involves  two  rotations;  the  first  about  the  Zj-axis  through  the  horizontal  flight  path 
angle,  o,  and  the  second  about  the  intermediate  y-axis  through  the  vertical  flight  path 


angle  y  (Figure  13.6).  The  rotation  matrix  required  to  transform  a  vector  X  expressed 
in  inertial  coordinates  into  the  intermediate  system  coordinates  can  be  written  as: 


R,{o} 


cosa  sina  0 
-sine  COSO  0 
0  0  1 


And  the  rotation  matrix  from  the  intermediate  system  to  the  flight  path: 


^2(Y) 


cosy  0  -siny 
0  10 
siny  0  cosy 


Multiplying  the  two  matrices  results  in  the  desired  transformation  matrix  to  the  flight 
path  coordinate  system  from  the  inertial: 


TFPT  =  R^iy)  R^ia) 


cos  y  COSO 
-sino 
_siny  COSO 


cos  y  sino 

COSO 

siny  sino 


-siny 

0 

cos  y 


(13.7a) 


Transposing  this  matrix  results  in  the  transformation  to  the  inertial  from  the  flight 
path  system: 


siny  COSO 
siny  sino 
cos  y 


cos  y  COSO 
-piFP  _  cosy  sino 

-siny 


-sino 
cos  o 
0 


(13.7b) 
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13.4.7  INERTIAL  FROM  INERTIAL  NAVIGATION  UNIT 
13.4.7.1  C-141  or  F-16  AIRCRAFT 

For  the  TPS  C-141  aad  F-16  aircraft,  the  transformation  to  the  inertial  reference 
frame  from  the  navigation  reference  frame  includes  two  rotations  (References  2  and 
3);  the  first  about  the  axis  through  an  angle  of  180  degrees  to  align  the  z-axes 
(note  that  Zj^u  is  up),  and  the  second  about  the  Zj^u-axis  through  the  wander  angle, 
ttiNu,  as  shown  in  Figure  13.11.  The  rotation  matrix  required  to  transform  a  vector^ 
expressed  in  navigation  system  coordinates  into  the  intermediate  system  coordinates 
can  be  written  as: 


R,{180°)  = 


10  0 
0-10 
0  0-1 


And  the  rotation  matrix  from  the  intermediate  system  to  the  inertial: 


R.  (a. 


cosa^^ 

-sina^^ 


sina^ 


cosa 


INU 


0 

0 

1 


Figure  13.11  C-141/F-16  Inertial  Navigation  Unit  to  Inertial  Coordinate 
System  Transformation  through  the  Wander  Angle,  Uinu 


Multiplying  the  two  matrices  results  in  the  desired  transformation  matrix  to  the 
inertial  system  from  the  navigation  system: 
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Ti,iM7  =  R^(a,^)  R,{180°) 

cosa^^  -sina^^  0 

=  -sina^^  -cosa^^  0 

0  0-1 


(13.8a) 


13.4.7.2  F-15  Aircraft 

For  the  F-15  aircraft,  the  transformation  to  the  inertial  reference  frame  from  the 
navigation  frame  is: 


yi,  INU 


10  0 
0  10 

0  0-1 


(13.8b) 


This  matrix  simply  provides  the  multiphcation  of  the  Zj^u  component  of  the  given 
vector  by  -1  (Reference  4),  converting  the  navigation  left-handed  system  (North,  East, 
Up)  to  a  right-handed  inertial  system  (North,  East,  Down). 


13.4.8  BODY-FIXED  FROM  THRUST 

The  transformation  from  the  body-fixed  to  the  thrust  reference  frame  includes  a  single 
rotation  about  the  Zg-axis  through  the  angle  Ot  followed  by  a  second  rotation  about  the 
intermediate  y-axis  through  an  angle  Yt,  as  shown  previously  in  Figure  13.7.  This 
results  in  a  transformation  matrix  similar  to  that  found  for  the  inertial  to  flight  path 
transformation: 


'pi’B  _ 


COS  YrCOSOy 
-sinGy 
sin  Yy  cosGj, 


cos  YySinGj, 

cos  Gy 

sin  YySinGy 


-sinYr 

0 

cos  Yr 


(13.9a) 


But  what  we  really  need  is  the  inverse  transformation  to  go  to  body-fixed  from  thrust 
coordinates: 


'psr  = 


COS  YrCOSGy 
COS  YySinGy 

-sin  Yt 


-sinGy  sinYrCosGy 
cos  Gy  sinYrSinGy 
0  cos  Yr 


(13.9b) 


With  the  gross  thrust  vector  aligned  with  the  Xy-axis,  this  transformation  results  in: 


cos  Yr  cos  Gy  ' 

II 

0 

^  ^ 

COSYr  sinGy  >. 

B 

^  0 

-Fg  sinYr 

(13.10) 
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13.5  GENERAL  VELOCITY  RELATIONS 

13.5.1  LINEAR  VELOCITY 

The  true  airspeed  of  an  aircraft,  V^,  calculated  from  the  air  data  system  represents  the 
magnitude  of  the  aircraft  flight  path  vector  relative  to  the  surrounding  air  mass.  This 
flight  path  vector  is  directed  along  the  wind  coordinate  system  x-axis  as  shown  iu 
Figure  13.5.  As  long  as  the  air  mass  is  not  moving,  the  true  airspeed  can  be  related 
directly  to  the  ground  speed,  Vg,  through  the  transformations  defined  in  section  13.4.5. 
However,  if  the  air  mass  is  moving,  the  calculation  of  ground  speed  must  include  the 
velocity  of  the  air  mass  relative  to  the  ground,  V^.  Therefore,  the  general  relationship 
between  the  aircraft  velocity  in  the  inertial  frame  and  the  aircraft  velocity  relative  to 
wind  can  be  written  as: 


' 

0 

^  ^ 

I 

0  , 

w 

kJ 

[  c0  c'F  ca  cP  +  { s 0  cT  -  c3> sT)  s p 
+  (c<E)s0  c'F  +  sa  cPJ  + 

[c0  s'Pca  cP  +  (s<I>s0  3^?+  c<I>c'P)  sP 

+  (c^s0  sT-s<I)c'P)  sacp]  +  v„ 

Vt  [-S0  ca  cp  +  sOc0  sp  +  c^c0  sa  cp]  +v„ 


If  an  aircraft  is  equipped  with  an  INU,  the  inertial  velocity  can  be  measured  directly, 
allowing  the  wind  velocity  components  to  be  calculated  from  the  above  equations. 
Notice  that  the  ground  speed  is  simply  the  velocity  in  the  inertial  horizontal  plane, 
and  is  given  by: 

=  \/V  ^  (13.12) 

Additionally,  the  aircraft  actual  climb  rate,  ti,  is  simply  the  negative  of  the  vertical 
component  of  velocity 

*  =  (13.13) 

Therefore,  from  equation  13.11, 

h  =  v^[  {s0  ca  -  c$c0  sa)  cp  -  sOc0  sp]  - 


(13.14a) 
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Measurements  of  climb  rate  from  an  aircraft  air  data  system,  H^,  must  be  corrected 
since  they  differ  from  actual  (geometric)  climb  rates  due  to  deviations  in  ambient 
temperature  from  standard  day  values.  Thus,  equation  13.14a  can  be  rewritten  as: 

h  =  =  V,  [  (s0  ca  -  c<I)c0  sa)  cP  -  s^c0  sPl  - 

“  ‘  K  K  ...  (13.14b) 

“  ®  YinertiaJ ) 

where  is  the  actual  ambient  temperature  at  the  test  altitude,  is  the  standard 
day  ambient  temperature  corresponding  to  the  test  altitude,  YAirwass  aircraft 

flight  path  vector  relative  to  the  air  mass,  and  Ymerciai  aircraft  flight  path 

vector  relative  to  the  inertial  (flat  Earth)  frame. 


13.5.2  ANGULAR  VELOCITY 

Based  on  the  angular  relationships  among  the  aircraft  coordinate  systems  defined 
previously,  we  can  now  define  the  angular  velocity  vector  between  any  two  systems. 
First  consider  the  angular  velocity  between  the  inertial  and  body-fixed  coordinate 
system.  Three  rotations  occur  through  the  Euler  angles,  and  the  angular  rates 
associated  with  each  rotation  can  be  added  vectorially.  Thus,  the  angular  velocity  of 
the  body-fixed  axes  with  respect  to  the  inertial  axes  is: 


=  'F  +  0  J2  +  6 

Expressing  this  vector  in  terms  of  body-fixed  coordinates: 


(13.15a) 


or 


{0)}i 


0 

0 

6 

0 

^  (O)  “ 

0 

^  ^ 

0  ► 

T 

I 

0 

2 

0 

{COjj 


-'Fsin0 

'  > 

0 

f  •  ^ 

<I) 

'Psin3)cos0 

^  ^ 

0  cos^ 

>  +  . 

0  ^ 

Ycos<I>cos0 

B 

-0  sinO 

B 

0 

(13.15b) 
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But  the  components  of  this  angular  velocity  vector  expressed  in  terms  of  the  body- 
fixed  coordinates  are  defined  as  the  inertial  angular  rates  P,  Q,  and  R,  where: 


0)„  s  Inertial  roll  rate,  P 
=  Inertial  pitch  rate,  Q 
=  Inertial  yaw  rate,  R 

(13.16a) 
(13.16b) 
(13.16c) 


Thus,  the  inertial  body  system  angular  rates  are  given  by: 

P  =  6  -  'PsinG 

Q  =  vp  sin<S)  COS0  +  0cos<l> 

R  = 'P cosO  cos  0  -  0  sin^> 


Note  that  if  we  wish  to  instead  relate  the  angular  velocity  of  the  inertial  frame  with 
respect  to  the  body,  we  need  only  multiply  the  previous  results  by  (-1),  since 

(13.17) 

Now  consider  the  angular  rates  between  the  body-fixed  and  wind  coordinate  systems. 
The  equation  for  the  angular  velocity  of  the  wind  axes  with  respect  to  the  body-fixed 
axes  can  be  expressed  as 


d®" =  -aj. 


(13.18a) 


Notice  the  negative  sign  of  a  is  due  to  the  definition  of  a;  a  positive  rotation  from  the 
stability  to  the  body  axes.  Our  transformation  took  us  the  opposite  way.  Now, 
expressing  equation  13.18a  in  terms  of  body  coordinates; 


or 


{0)}b 


BW 


"  'I 

''  ^ 

0 

0 

-a 

>  +  ’J'BS  ^ 

0 

0 

B 

P 

{coir  =  < 


-psina 

-a 

pcosa 


(13.18b) 
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Again,  the  angular  velocity  of  the  body-fixed  system  with  respect  to  the  wind  system 
is  simply  the  negative  of  the  above  results 

(13.19) 

Note  that  we  can  develop  any  angular  velocity  vector  between  two  reference  frames 
and  express  it  in  terms  of  a  desired  coordinate  system  based  on  similar  analysis. 

Recall  that  to  apply  Newton’s  2“'*  law,  time  derivatives  can  only  be  taken  in  an 
"inertial"  reference  frame.  The  primary  significance  of  the  angiilar  velocity  vector  is 
that  it  allows  the  calculation  of  an  inertial  time  derivative  of  a  vector  expressed  in  a 
rotating  frame  to  be  performed  in  the  rotating  frame,  as  illustrated  in  the  following 
relation: 


^  dAji  _  ^ 

dt  dt 


+  0) 


NR,R 

R 


xA^ 


(13.20) 


where  Ajj  is  an  arbitrary  vector  expressed  in  the  rotating  frame.  For  example, 
consider  the  true  velocity  vector,  which  is  defined  in  the  wind  reference  frame.  The 
time  rate  of  change  of  this  vector  in  the  inertial  reference  frame  can  be  expressed  as: 


J  _  W  _ 
d  V.  d  V. 


dt 


dt 


X 


Q 7  =  } 


where 
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13.6  GENERAL  FORCE  RELATIONS 

13.6.1  AIRCRAFT  EXTERNAL  FORCES 

The  external  forces  acting  on  an  aircraft  include  the  gravitational  force,  aerodynamic 
forces,  and  propulsion  system  forces.  Each  of  these  forces  is  defined  in  one  of  the 
coordinate  systems  presented  in  section  13.3,  as  shown  in  Figure  13.12. 


L 


w 


Figure  13.12  Aircraft  External  Force  Definition 

The  gravitational  force  vector  is  along  the  Zj-axis,  and  is  given  by: 

'^grav  =  W  =  lug  =  mgk^  (13.21) 

where 

m  =  aircraft  instantaneous  mass 
g  s  gravitational  acceleration, 

assumed  constant  (=  32.2  ft/s^  or  9.81  m/s^) 
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The  aerodynamic  force  vector  is  defined  in  the  stabiUty  coordinate  system  as  follows: 


where 


F. 


=  D  +  Y,  +  L  =  -D 1+  Y. 


J, 


-  Lk^ 


D  =  Aerodynamic  Drag 

=  Aerodynamic  Side  Force 
L  =  Aerodynamic  Lift 


(13.22) 


The  propulsion  system  forces  include  the  gross  thrust  and  engine  drag.  The  gross 
thrust,  Fg,  is  defined  in  the  thrust  coordinate  system  while  the  engine  drag,  F^,  is 
defined  in  the  stability  system.  Thus,  the  propulsion  system  force  vector  can  be 
written  as 

^p.op  =  Fjr  -  fJs  (13.23) 


13.6.2  GENERAL  FORCE  EQUATION  OF  MOTION 

From  Newton’s  2nd  Law,  the  sum  of  the  external  forces  is  equal  to  the  time  rate  of 
change  in  momentum: 

_  J  ,  I’w  _ 

Ef,,,  =  ^  imv)  (13.24a) 


Even  though  the  vehicle  mass  is  constantly  changing,  we  can  consider  it  as  a  constant 
in  the  equation  for  an  instant  in  time.  This  results  in 


Ef 


ext 


=  a 


(13.24b) 


Inserting  all  external  forces  acting  on  the  aircraft  yields 


0 

0 

'  + 

>  +  ^ 

0  ’ 

mg 

I 

-L 

3 

.0  ^ 

(13.25) 


Expressing  this  equation  in  terms  of  the  stability  coordinate  system,  assuming  the 
inertial  accelerations  are  originally  expressed  in  the  inertial  frame  (possibly  measured 
from  an  INU)  and  do  not  include  gravitational  acceleration,  produces: 


V 
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m  1 


0 

0 

[mg 


-D-F^ 

-L 


+  T^T  \ 


'-mg  ( s0  ca  -  cO  c6  sa)  -  D  -  { sy^  sa  -  cy^  cg^  ca ) 

=  <  mg  sO  c0  +  +  Fg  cYr  sa^ 

mg  {s0  sa  +  c4>  c0  ca  )  -  L  -  Fg  { sy^  ca  +  cy^  co^  sa)  ^ 

(13.26a) 

Equation  13.26a  represents  the  most  general  form  of  the  force  equation  expressed  in 
stability  coordinates.  If  we  assume  that  the  thrust  vector  is  fixed  at  an  angle  i^  with 
the  x-body  in  the  x-z  body  plane,  then  this  equation  simplifies  to 


-mg  ( s6  ca  -  c4)  c6  sa)  -  D  -  F^  +  Fg  c  [a  +  i^)  ' 
mg  sO  c0  +  , 

mg  (s6  sa  +  cO  c0  ca  )  -  L  -  Fg  s  {a  +  i^) 


(13.26b) 

If  we  consider  the  following  additional  assumptions  for  an  aircraft  in  steady,  wings 
level  flight: 


1.  Zero  Roll  (O  =  0) 

2.  Thrust  vector  fixed  along  x-body  axis  (Yt  =  Ot  =  0) 

3.  Zero  winds 

4.  Zero  inertial  accelerations  (a,  =  ay  =  a^  =  0) 


equation  13-26b  reduces  to 

(-mg  sin  (0  -  a) 


0  = 


-  {F^  +  D)  +  Fg  cosa' 


mg  cos  ( 0  -  a)  -  L  -  Fg  sina 


S 
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and  the  classical  equations  emerge: 


Fg  cosa  -  -  mg  siny  =  F^  -  W  siny  =  D  (13.27a) 

=  0  (13.27b) 

mg  cosy  =  W  cosy  =  L  +  F^  sina  (13.27c) 

where 

=  net  thrust 
W  =  aircraft  weight 


The  accelerations  imposed  on  an  aircraft  can  also  be  measured  with  body 
accelerometers.  However,  these  measured  accelerations  include  the  gravitational 
acceleration  acting  on  the  aircraft.  Therefore,  equation  13.25,  expressed  in  the 
stability  reference  frame,  must  be  rewritten  as: 


m 


a 


accel 


) 


-d-f; 

,  H-  JST  < 

0  > 

-L 

S 

.0  ^ 

D  -  F^  -  Fg  ( syj.  sa  -  cy^  ca^  ca) 

+  Fg  cy^  SQy 

-  L  -  Fg  {sy^  ca  +  cy^  cc^  sa) 


(13.28a) 


Again,  if  we  assume  that  the  thrust  vector  is  fixed  at  an  angle  i^  with  the  x-body  in 
the  x-z  body  plane,  then  equation  13.28a  simplifies  to 


m 


-  D  -  F^  +  Fg  c  {a  +  i^) 
Ya 

-L  -  Fg  s  {a  +  i^) 


(13.28b) 

Equations  13.26  and  13.27  are  valid  for  application  to  any  stable,  quasi-steady,  or 
dynamic  performance  test  method.  In  the  second  half  of  this  course,  we  will  use  the 
same  principles  presented  to  complete  the  development  of  the  aircraft  moment 
equations  of  motion  for  application  to  flying  qualities  testing. 
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